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1. INTRODUCTION 

Increased capacity demands, improved throughput, and additional use cases for wireless access 
where all things will be connected in a seamless fashion have steered a new vision of mobile broadband 
communications [1]-[4]. The 5th generation (5G) entails new technologies, which have to fulfill many 
requirements such as the high data rates per user, uniformly good service throughout the coverage area and 
the ability to function under high mobility conditions, etc. The deployment of small cells and exploitation of 
higher frequencies are known as 5G key technologies, nevertheless they are unable alone to address these 
needs. 

Massive MIMO (Multi Input and Multi Output) is the new and the most promising direction in 
mobile access that paves the way for 5G [5]. This emerging technique is an eventual and worthwhile form of 
the Multi User MIMO called also Large scale antenna system, very large MIMO, or hyper MIMO [6]-[14]. It 
scales up conventional MIMO by possibly of degrees of magnitude compared to the current state-of-the-art. 
The Massive MIMO technique could constitute a network of cells, or it could be an isolated (single-cell) site 
where antenna arrays with a very large number of physically small, non-directive antennas are employed to 
serve a multiplicity of autonomous users (terminals with a single antenna) in the same time-frequency 
resources. 

Generally, Massive MIMO Systems operate in Time Division Duplex (TDD) mode and rely on 
reciprocity between the uplink and downlink channels. It is an obvious choice for two good reasons. Firstly, 
the information about the channels is required only by the BS to process the antennas coherently. Secondly, 
the uplink estimation overhead is proportional to the number of terminals but independent of service antennas 
number. However, the uplink resources required in the Frequency Division Duplex mode (FDD) to inform 
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the BS about the channel responses would be drastically bigger and proportional to the number of antennas at 
the base station [15]. In spite of that fact, the FDD operation may be possible in certain cases [16]. 

On TDD massive MIMO systems, an orthogonal uplink pilot sequence is assigned to every terminal. 
In fact, the limited duration of the coherence interval is the main reason that the orthogonal pilot sequences 
are normally reused by users in neighboring cells. In [17], for a typical operating scenario, the maximum 
number of orthogonal pilot sequences, in a one-millisecond coherence interval, is estimated to be about 200. 
The accessible supply of orthogonal pilot sequences in a multi-cellular system is typically expended. Hence, 
the necessity of reusing non-fully orthogonal pilot sequences across all the cells. Unfortunately, the multi-cell 
massive MIMO system would be plagued by the so-called pilot contamination. In particular, a linear 
combination of terminals’ channels having the same pilot sequence contaminates the estimated channel 
obtained in the service array. 

The pilot contamination has a detrimental impact on the actual achievable spectral and energy 
efficiencies in real systems [18]-[20]. As a result, considerable research efforts have been spent in the last 
years towards alleviating pilot interference in massive MIMO networks. It has been demonstrated in [21] that 
the effect of pilot contamination can be exterminated using a Bayesian channel estimation ensuring certain 
conditions of the channel covariance. The main idea in [21] is to assign the same pilot to users revealing 
maximum orthogonality of signal subspaces (with the same Direction Of Arrival DOA). Besides, all 
covariance matrices of all channels must be learned previously which incurs too much coordination. Another 
nonlinear channel estimation technique is proposed in [22] as a blind pilot decontamination method. It 
consists of exploiting the difference between the channel gains of the desired users and the interfering users 
in order to discriminate users in the amplitude domain. This approach was considered in [23] for a more 
sensible physical channel model. The authors discussed the ability of subspace method to identify the desired 
and interfering users. It was proven that the subspace approach works well within the constraint that the 
interference channel is received with a power level sufficiently lower than that of the desired channel. 

In effect, most of the previous works considered the Elevation parameter and/or the constraint of a 
larger power ratio between the intended users and the interfering ones. Nonetheless, in real channels, the 
Elevation parameter (1D DOA) is a finite number and can be identic in some scenarios such a rich scattering 
environment. As well, the constraint of power ratio is a condition hard to guarantee for some edge-of-cell 
users as they may receive their signals with the same Elevation angle (El) and with the same power. 

Therefore, we aim in this work to propose an effective solution to the pilot contamination problem 
for all cell users and more precisely taking into account those edge-of-cell users. Our idea behind is to state 
and resolve the pilot contamination problem as a source separation problem using time, frequency, and space 
domains. The mixing matrix compromises both the interfering and desired sources propagation channels. The 
latter is defined by the sources spatial signatures bearing angle that will be recovered by a DOA estimation 
technique. Assuming that 1D DOA does not allow separating the intended users and the interfering ones 
since the signals of two users may be received with the same direction, we propose to include both El and Az 
angles to separate the users of each cell. We assume that the probability of two users having the same values 
for both elevation and azimuth angles is very small. 

In this paper, we exploit the concept of STFD [24]-[26]. It has been used to provide high resolution 
and enable separations of signals corrupted by interference occupying the same frequency band and/or the 
same time segment. It offers an improving signal selectivity over approaches using the covariance matrix. 
The correlation matrix is a low-dimension representation of STFD, the reason that STFD-based methods 
offer better performance than conventional covariance matrix based approaches. Afterward, we apply a 
nonunitary Joint diagonalization (JD) method to a combined set of Spatio-temporal correlation matrices, 
which is selected using a Time-Frequency TF point detector cited in [26] to provide an improved estimate of 
the signal and noise subspaces over the standard Singular Value Decomposition (SVD) technique. Finally, 
the Multiple Signal Classification (MUSIC) technique is used to determine the values of El and Az angles. 
The result is an improved sources separation using the parameters: time, frequency, and space variables. 

The outline of this paper is as follows. We describe in Section 2 the problem formulation and derive 
the concept of the STFD representations for a Multi-Cell Massive MIMO system. We discuss the challenges 
for selecting TF points and present the method applied for selecting the TF points. We introduce also the NU- 
JD approach used for estimated the noise subspace that is used in its turn to determine the azimuth and 
elevation angles using MUSIC spectrum. Simulations results and conclusion are presented in Sections 3 and 
Section 4. 
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2. PROBLEM FORMULATION 
2.1. System model 

We will consider the uplink in multi-cell multi-user MIMO communication system with L cells as 
shown in Figure 1, with full spectrum reuse. There are K single-antenna users in each cell simultaneously 
served by their Base Station (BS) equipped with a Uniform Rectangular Array (URA) as shown in Figure 2. 
The cellular network operates in time-division duplexing (TDD) mode and the channel is assumed to be a 
block-fading channel, which remains constant over a coherence time of C symbol periods, and changes 
independently from one coherence time to another. The same pilot pool is reused among all cells and the 
pilot sequences are assumed orthogonal to each other in order to avoid intra-cell interference. 
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Figure 1. The concept of pilot contamination in a multi cell massive MIMO system 
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Figure 2. System model of M x N URA configuration 


In this paper, we consider that the carrier frequency is 3GHz and the URA is selected with the 
consideration of limit space. The antenna array at the base station is placed in the X- Y plane with M antenna 
elements in each column and N antenna elements in each row. The spacing between adjacent antenna 
elements is assumed to be half wavelength (d < A/2) in both two dimensions. The DOA is depicted by 
elevation and azimuth angles, given by 8 € {0°,90°} and € {—90,° 90°} respectively [27], [28]. 

The so-called steering vector quantifies the array response in a given direction. Hence, the channel 
from the j user in the i'"cell to the intended BS is denoted by hji. We adopt the Correlation Channel Model 
[29] given by: 


hj; = q%iaji (8j Hii) (1) 
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a. aji: is the steering vector corresponding to the angles of arrival, the elevation 9j; and the azimuth 
Qþjiassociated with the j user. 
b. aji~CN (0, Bi) : is the channel gain. 


Moreover, it is assumed that all 6;;, bj; et aj; are independent over user index j and cell index i. The 
MN x 1 manifold aj (9j;, ji) is given by 


2nd amd T 
= —— sin(6;;) ——(N-1) sin(6;;) 
ayi (8 ji, ji) = vec{|1,e x ye A it 


amd- 2nd J 
@[1, ex Ni) ge MD sin(bji)}T} (2) 


Here after ® et vec{. } operators the Kronecker product and the vectorized operation respectively. 
Let H; 2 (Aiii ..., Axjhy;) be the MN x K effective fast-fading channel from K users in the it” cell to the 
intended BS, where: 


a. Ajj is the steering vector of the j th user in the i” cell; 

b. hy € C*** consists of iid CN (0,1). 

The MN x C signal received by the 1-th BS during C consecutive symbol intervals may be written 
as: 


Y =H,X, + XL, HiX; +B (3) 


where B is an additive noise vector whose entries are considered generally as stationary, temporally and 
spatially white, zero mean random processes, and independent of the sources signals. 
The Equation (3) can be simplified to the following equation: 


Y=HX+Z (4) 


X=X,e CK*C is the transmitted signals from the intended cell sources through the channel 
H =H, € CMN*K and the impairment process Z € CMN*C is decomposed into white noise B and interference 
from L- 1 neighboring cells 


where the data X; € C-YK*C is transmitted from interfering cells and received in the cell of interest 


through the channel H; € CMN*@-DK | 
Now, we consider the nonlinear model of the studied system defined in Equations (4) and (5), the 
STFD of the dimension MN x MN is denoted as 


Dyy(t, v) = HDxx(t, v)H" + HDyy, (t, v)H + HyDx,x(t, v)H™ + HyDy,x, (t, V)HF + HDxg(t, v) + 
Dpx(t, v)H# + H,Dx p(t, v) + Dex; (t, v)HE + Dpp(t, v) (6) 


where: 


a. Dxx(t,v) represents the STFD of the intended cell sources, 

b. Dx,x,(t v)is the STFD of the interference 

c. Dyx,(t,v) and Dy,x(t, v)are the crossed STFDs between the sources and interference from L — 1 
neighboring cells 

d. Dyg(t,v) and Dpgx(t, v) symbolize the crossed STFDs between the intended cell sources and 
noise, 

e. Dx p(tv) and Dgx,(t,v) are the crossed STFDs between the interference from L—1 
neighboring cells and noise. 

f. Dyn(t, V) represents the STFD of the Noise. 
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Under the assumption of white noise, the expectation of the STFD matrices between the intended 
sources signals and the noise Dxg(t, v), and Dgx(t, v) and between the interference from L neighboring cells 
and noise Dx,g(t, v) and Dgx,(t, v) vanish, i.e. 


E[Dxp(t, v)] = 0, E[Dgx(t, v)] = 0, 


E[Dx,8(t,v)| = 0,E[ Dgx,(tv)] = 0 
The interference is not independent of the data X, hence the STFD of the observations becomes: 


Dyy(t,v) = HDyx(t,v)H" + HDy,,(t,v)H} + H,Dy,x(t,v)H4 + H,Dy,x,(t, v) Hf + 
Dgg(t, v) (7) 


2.2. Pilot decontamination using time-frequency El-Az estimation 
The proposed method for pilot contamination problem can be summarized in Figure 3 follows 


Figure 3. The Organogram of the proposed method NU-TF-MUSIC El-Az estimation 


El-Az estimation can only be materialized if appropriate TF points are selected in the STFD 
matrices formulation as shown in Equation (7). 


2.2.1. Time-frequency points selection 

We use the same TF points detector proposed for the instant context. It concerns finding useful TF 
points through inspecting the time frequency points where the STFD observations matrices are rank one. For 
this reason, a Singular Value Decomposition (SVD) is used leading to 


Dyy(t,v) = (U(t,v) A(t, v) V(t v)") (8) 
with V(t, v) et U(t,v) are two MN x MN unitary matrices, and A(t, v) = diag{5(t, v)} is a diagonal matrix 
whose elements are positive. Let 6(t,v) = (5,(t,v), ...,5,(t,v))” and nd assume that the singular values are 


arranged in a decreasing order 6,(t,v) > 62(t,v) > = > 5x(t,v) = 0. Therefore, we choose a point (t, v) 
and check if a matrix is rank one according to the following two conditions. 


§,(t,v) > &, 


M8 v) < £z, (9) 


where €,(resp. €2) is a positive constant sufficiently large (resp. small). The set of Ny(with Ny E N*) of 
matrices M which is being built, all admit a particular structure such that 


Dyy(t, v) = H Dxx(t, v)H" + HDxx, (t, v)HË + HyDx,x(t, v)H" + HyDx,x, (tv) Hf (10) 


with Dxx(t, v),Dxx,(t,v), Dx,x(t, v) and Dx,x,(t, v)are diagonal matrices with a single non-zero term on the 
diagonal. 


2.2.2. Joint diagonalization 
We consider the set M built using the detector that has just been presented in the previous section. 


M = {M; = Dyy(ti, vi) E CMxM yi = 1, wy Nm} 
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All matrices of the set M admit the following decomposition: 


M; = ED,E® or D; = FM;F"vi € {1, ..., Nm} (11) 
where: 
di 0 0 
D; = i = l vi {1, .., Nu} 
0 = 0 dikk 


Are diagonal matrices of dimension LK x LK and dij i€ {1,..,Nm}j=ke {1,...,LK} are the auto- 
terms. The matrix E is a full rank matrix of dimension MN x LK and F is its pseudo-inverse matrix. 

The problem of nonunitary joint diagonalization is to estimate the separation matrix A only from the 
matrices of the set M. We consider the following criterion: 


Cp (â) = DMM lOffdiag{AM,A"} ||" (12) 


where ||. ||p denotes the Frobenius norm and the Offdiag{.}operator extracts the diagonal elements of a 
matrix to build a zero diagonal matrix. In order to find the estimated A, an estimate of the steering matrix A, 
we need to minimize the sum of squared norms non-diagonal elements of AM,A". This is done as in [29] 
through an algebraic solution based on computing the exact Hessian complex matrices. 


2.2.3. El-Az estimation 

The El-Az estimation is accomplished using NU-TF-MUSIC, which looks for peaks in a spatial 
spectrum computed for every pair of angle in a region of interest. To describe NU-TF-MUSIC, we denote Ê, 
as the noise subspace of the estimated steering matrix obtained by jointly diagonalizing the STFD matrix. 
The estimation technique MUSIC estimates jointly the El and Az angles by determining the LK values of 0 
and @ for which the following spatial spectrum is maximized, 


1 


TF = 
fuusic 8, p) E [a(0,o)HE,EHa(0,o)]7 (13) 


Once the noise subspace is obtained, the El and Az angles are determined by locating the LK peaks 
of the spatial spectrum. After all, the pilot contamination problem is cast as sources separation problem as we 
obtained the El and Az angles of every distinguished user. The BS can establish its users’ channels and edge- 
of-cell users’ channels without falling back on the constraint of power control to ensure that the intended 
users’ power is upper than the interfering users power. Consequently, the effect of pilot contamination is 
eliminated. 


3. PERFORMANCE EVALUATION 

In this section, we evaluate the performance of the proposed approach in a multi-Cell Massive MIMO 
context by conducting Monte Carlo simulations. In the simulation, we consider two hexagonally shaped 
adjacent cells in the network. Each base station has a total number of 128 antenna elements, which is 
considered as one possible dimension of massive MIMO [30]. The 128 antenna elements are extended over a 
M x N dimension forming a URA configuration, with half wavelength antenna spacing. There are four users 
per cell using the same pilot sequences and having the same elevation angle values where 
Bce = (5°, 15°, 25°, 35°} Ocens = (5°, 15°, 25°, 35°}, Del = {0°, 8°, 20°, 90°} and 
Dcell2 = {10°, 30°, 40°, 60°}. We consider the root-mean-squared error (RMSE) as a performance metric to 
measure the effectiveness of El-Az estimation and sources separation, which is defined as: 


RMSE = [DES Ene (x - Ox) + (Ox - ox) | (14) 


where 6), and ¢;, are the estimation of elevation angle ®, and the estimation of azimuth angle $, respectively 
at the n Monte Carlo trial. LK is the estimated source number and N,,. represents Monte Carlo simulation 


Int J Elec & Comp Eng, Vol. 8, No. 5, October 2018 : 3609 — 3619 


Int J Elec & Comp Eng ISSN: 2088-8708 O 3615 


times. We first show the spectrum of the proposed method, NU-TF-MUSIC, based only on the elevation 
angle for eight sources, where the four interfering sources and the four intended sources have the same 
values. The simulation is conducted with a coherence interval of C = 128, and a value of SNR = 15dB. 
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Figure 1. NU-TF-MUSIC spectrum for a SNR value of 15dB 


As we can observe in Figure 1, the spectrum exhibits just four peaks. One peak for each elevation 
angle({5°, 15°, 25°, 35°}). We cannot decide which user location is displayed in the figure. We can conclude 
that 1D DOA (only the elevation) does not allow discriminating users of each cell. In the following, we plot 
the NU-TF-MUSIC spectrum based on the elevation and azimuth angles of the proposed approach given by 
Error! Reference source not found.) for eight sources, where four interfering sources have the same 
elevation angle values as the four intended sources. 
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Figure 5. NU-TF-MUSIC spectrum for a SNR value of 15dB 


Figure 5 depicts the angles estimation performance of the proposed algorithm with the 16 x 8 array 
configuration, a coherence interval of C = 128, and a value of SNR = 15dB. The results show that 
elevation and azimuth angles can be clearly observed. The spectrum has been computed for all values of 0 
ranging from 0° to 90° and ¢ ranging from —90° to 90°, with an angular step size of 1°. NU-TF-MUSIC 
spectrum exhibits a peak at each source location (5°,0°), (5°, 10°), (15°, 8°), (15°, 30°), (25°, 20°), 
(25°, 40°), (35°, 90°) and (35°, 60°). 

Afterward, we present the RMSE of the proposed approach versus different SNR values and we 
compare the result for two types of the array configuration. 
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Figure 6. Comparison of RMSE of NU-TF-MUSIC and Conv-MUSIC versus SNR for different types of 
array configurations 


We show in Figure 6, the comparison of RMSE of NonUnitary Time Frequency MUSIC method, 
denoted by NU-TF-MUSIC, and Conventional MUSIC method, denoted by Conv-MUSIC, versus different 
SNR values for the four interfering sources having the same elevation angle values as the four intended 
sources. The value of RMSE of both methods and for different types of arrays is under | with SNR ranges 
from OdB to 30dB. As it has been pointed out as the SNR increases, the performance of the two methods 
ameliorates. The performance of each array configuration for the two methods is illustrated by a curve as 
mentioned in Figure 6. The 16 x 8 array configuration has acted in excellent quality for both NU-TF-MUSIC 
and Conv-MUSIC. Comparing the two array configurations, it is interesting to notice that the number of 
antenna element along X and/or Y axes contributes to improving the accuracy of angles estimation. The 
reason for which we observe a considerable difference between the RMSE values ranging from 0 to 15dB for 
the two configurations. This open issue will be discussed with researchers of our Laboratory working on 
Antenna subjects. 

Next, we display also the performance of the proposed method with another array geometry to see 
the impact of antenna geometry on our method. We choose the Uniform Circular Array UCA represented by 
the steering vector at any direction as in [31], [32] 
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Figure 7. NU-TF-MUSIC spectrum for a UCA configuration and for a SNR value of 15dB 
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Figure 7 illustrates the angles estimation performance of the proposed algorithm for a UCA 
configuration, a coherence interval of C = 128, and a value of SNR = 15dB. The elevation and azimuth 
angles can be clearly observed also for UCA configuration. We can conclude that our proposed method 
performs well for the URA and UCA configurations. 

In the following, we choose the 16 x 8 array configuration to compute the rest of simulations. We 
evaluate the RMSE of NU-TF-MUSIC and Conv-MUSIC versus different SNR values for two number of 
samples values. 
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Figure 8. Comparison of RMSE of NU-TF-MUSIC and Conv-MUSIC versus SNR for different values of 
coherence time C 


Figure 8 presents the NU-TF-MUSIC and Conv-MUSIC with an array configuration of 16 x 8, and 
two different values of coherence interval C. Obviously, as we get more samples to estimate the steering 
matrix A, the angles estimation performance will enhance. All values of RMSE of the two methods are under 
0.9 with SNR ranges from 0dB to 30dB. The NU-TF-MUSIC method performs better than the Conv-MUSIC 
method for the two values of coherence interval C = 128 and 256. 

To get a more quantitative understanding of how El-Az estimation method performs in a multi-Cell 
massive MIMO system, we display another simulation to prove what a massive antenna number can bring an 
enhancement in our approach performances. 
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Figure 9. Comparison of RMSE of NU-TF-MUSIC and Conv-MUSIC versus the number of antenna 
elements for a SNR=15dB and for four intended sources and four interfering ones 
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We exhibit in Figure 9 the performance of the estimation error of the two methods versus the 
number of antenna elements (32, 64, 80, 96,112 and 128) covering the eight sources for a SNR value 
of 15dB. The values of RMSE is under 1.4 for the antenna elements number of 32, 64, 80 and 96 for both 
methods. We can emphasize that as the number of antenna elements increases (beyond 96), the performance 
of NU-TF-MUSIC and Conv-MUSIC get more efficient. The analysis and simulations in Figure 9 note that 
the Conv-MUSIC method outperforms our solution. However, to get ameliorated performance, it is necessary 
to raise the number of realisations of our method to attain a good match between experimental and theoretical 
results. In our case, it was difficult due to the limitation of our computing equipment. 

This study is a contribution to the ongoing discussions about the problem of pilot contamination 
outcoming from the reuse of pilots over cells in multi-cell Massive MIMO systems. We have addressed the 
problem of discriminating the target cell users and interfering ones especially those situated on edge of the 
cell. The originality of our solution lies in the fact that the pilot contamination is resolved as a source 
separation problem using time, frequency, and space variables. These results are in good agreement with the 
assumptions considered in the beginning stipulating that including 2D DOA (El and Az angles) would 
separate the users of each cell in contrast of other studies which considered a 1D DOA. One question still 
unanswered is whether this algorithm can be readily used in practice for real system with high dimensions. 


4. CONCLUSION 

In this paper, the pilot contamination problem, regarded as a bottleneck of Multi-Cell Massive 
MIMO Systems, was posed and solved as a sources separation problem over time, frequency and space 
domains. We have shown that including another parameter for DOA estimation (both elevation and azimuth) 
instead of 1D DOA (elevation only) brought more precision and succeed to discriminate users and enhance 
the sources separation in Multi-Cell Massive MIMO systems context. The idea fully exploits the concept of 
Spatial Time Frequency representations and their application to direction finding. We have introduced a new 
NU-TF-MUSIC method for El-Az estimation to distinguish between intended users and interfering ones. The 
result is an improved signal localization and separation using time, frequency, and space variables. 
Numerical experiments illustrated the accuracy and efficiency of the proposed technique for the Multi-Cell 
Massive MIMO systems. 
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